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Abstract. We study the problem of particles undergoing Brownian motion in an 
expanding sphere whose surface is an absorbing boundary for the particles. The 
problem is akin to that of the diffusion of impurities in a grain of polycrystalline 
material undergoing grain growth. We solve the time dependent diffusion equation for 
particles in a d-dimensional expanding sphere to obtain the particle density function 
(function of space and time). The radius of sphere expands as i? = (Kt)"^ where a is 
the growth exponent, K a constant and t the time. This allows the calculation of the 
survival rate or the total number of particles per unit volume as a function of time. 
We have obtained particular solutions exactly for the case where d = 3 and a parabolic 
expansion of the sphere. Asymptotic solutions for the particle density when the sphere 
expansion rate is small relative to particle diffusivity and vice versa are derived. 
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1. Introduction 

The kinetics of many physical systems can be described by first-passage properties of 
stochastic systems (Redner 2001). In particular, the problem of the survivability of a 
diffusing particle, or Brownian walker, in a confined domain ( "cage" ) with moving and 
absorbing boundary ("fence") has been a subject of several recent studies (Krapivsky 
& Redner 1996, Redner 2001, Bray & Smith 20076, Bray & Smith 2007a). 

Krapivsky & Redner (1996) evaluated the survival probability of a diffusing particle 
S{t) as a function of time t within a 1-dimensional expanding cage of length L{t) with 
absorbing boundaries by solving the standard diffusion equation (heat equation) in 
one dimension. Approximate methods were devised in the limit of slow (adiabatic 
approximation) and fast (free approximation) motion of the absorbing boundary. More 
specifically, Krapivsky & Redner (1996) considered L{t) ~ with a < 1/2 for adiabatic 
approximation and a > 1/2 for free approximation. Moreover, the case of a = 1/2, 
corresponding to marginally expanding cage, was theoretically assessed, see also Redner 
(2001). 

Bray & Smith (20076) calculated the exact asymptotic survival probability of a 1- 
dimensional Brownian particle, initially confined at a position x G {—L{t), L{t)) in the 
presence of two moving and absorbing boundaries by solving a Fokker-Planck equation 
for the system in the limit t oo. Subsequently, they (Bray & Smith 2007a) extended 
their method to treat the case of the d-dimensional system and in particular d = 3. 
Similarly, in the latter work, they considered the steady-state case (t oo) and assumed 
the absorbing sphere's size evolves as L{t) = Lq + vt, where Lq is the initial size and v 
is the constant velocity of the boundary. 

The aforementioned deliberations have relevance to the problem of diffusing 
particles (impurities) in a grain of polycrystalline material under grain growth (Forsberg 
& Massih 2007). Grain growth is a kinetic process at which the mean grain size of 
an aggregate of crystals increases during annealing at an elevated high temperature 
(Atkinson 1988). In this course, the larger grains are apt to increase in size at the 
expense of smaller grains which collapse, causing the total number of grains to decrease. 
Impurities in solid solution, prevalent in polycrystalline materials, which usually undergo 
diffusive motion at elevated temperatures, are in addition subjected to a moving or 
stretching medium due to grain growth. In particular, the process occurs in nuclear 
fuel (e.g., UO2), whereupon the fission product gases (e.g., Xe and Kr) diffusing in the 
grain can be subjected to a moving grain boundary. Grain boundary movement can 
sweep up fission gas atoms more rapidly than they could have arrived at the boundary 
by diffusion (Hargreaves & Collins 1976). This problem has been treated in (Forsberg 
& Massih 2007). Modelling grain growth phenomenon or evolution of cellular structure 
in general has been a subject of many recent investigations using various theoretical 
and numerical approaches, see review by Thompson (2001). For example, in (Niemiec 
et al. 1998, Gadomski et al. 2003) a different approach is used to describe the grain 
growth process and for which the exponents of the growth asymptotes are presented. 
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In this note, we solve the time-dependent diffusion equation for particles in a d- 
dimensional expanding spherical grain to obtain the particle density as a function of 
space and time. Equations of motion are treated in section [21 General solutions for 
particle density are provided in section [31 Special solutions are derived in section [H for 
the case of ci = 3 and a parabolic grain growth. An evaluation of particle survival is 
presented in section [Sj here we also include asymptotic solutions for particle density 
when grain growth rate is small relative to particle diffusivity and vice versa. We 
conclude the note in section [6l with a brief summary of the results. 



2. Equations of motion 



We consider the equivalent sphere model for material, meaning that the grains of 
the material are treated as a collection of spheres of uniform size characterized by a 
single equivalent radius, R = R{t), which is a function of time t. The particles (e.g. 
interstitials, impurities etc) migrate to grain boundaries by diffusion upon which are 
released from the system. The grain can contain traps (vacancies, pores, etc), which 
may capture the diffusing particle during its flight. The differential equation for the 
density of particles at position r in a (i-dimensional spherical grain at time t, C{r,t), is 
given by the law of isotropic diffusion, viz. 

dC{r,t) _ . 



dt 

subject to the conditions 
dC{0,t) 



DVlC{r,t), (1) 



Or 

cm),t) = o, (3) 

C(r,0) =Co, (4) 

where D is the effective diffusivity of particles in the grain, accounting for the presence 
of traps, = d'^/dr'^ + {{d — l)/r)d/dr, Co is the initial density of particles in the 
grain, and equation ([3]) defines the absorbing boundary (Dirichlet boundary condition). 
We should note that D, the effective diffusion coefficient of the particles in the grain, 
accounts for the effect of intragranular traps. As has been shown by Speight (1969), 
for stationary traps, it is related to the intrinsic diffusivity Dint of particle in the solid 
through D = pDint/{p + g), where p is the escape probability of the trapped particle 
and g is the corresponding capture probability, related to trap size, the number of traps 
n and the intrinsic diffusivity, e.g. for spherical traps, g = AirDintQ'n, where g is the trap 
mean radius (Ham 1958). Thus at any time, only a fraction of the particles is trapped, 
while the rest are able to diffuse out of the grain. 

The total number of particles per unit volume of d-dimensional spherical grain (V^) 
at time t, G(t), is given by 

G{t) = ^j\'-'C{r,t)dr. (5) 
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The total number of particles contained in the grain at time t is simply S{t) = Vd{t)G{t) 
with Vd{t) = n^^'^R'^{t)/T(d/2 + 1) being the volume of sphere at time t; while the total 
number of particles released from the grain is 

F{t) = V,{0)Co-Sit). (6) 

In general, we allow the particle diffusion and grain growth to occur simultaneously; 
hence we scale equation ([T]) to embed the time-dependent variable R{t) in the coordinate 
of the partial differential equation (space-time), namely 

dC{p,T) d'^C{p,T) id-1 Rr l9C(p,r) 



+ 



-F7P 



dr dp^ L p R i dp 

where we put 



(7) 



P^^. (8) 

ds 

with Rr = dR/dr and the conditions: 9C(0,r)/9p = 0, C(l,r) = 0, and C(p,0) = Cq. 
The total number of particles contained in the grain as a function of r is obtained by 
transforming equation ([5]) to 

g(T) = d [ x'^~^C{x,T)dx. (10) 
Jo 

Note that in our convention the variables are transformed to C{r,t) =>■ C{p,t), 
G{t) =^ (y(r) and i?(t(r)) = 7?.(r). Solving equation (I7j) for C{p,t), Q{t) can be 
determined through the evaluation of the integral in (fTOl) . 

Before discussing the solution of equation (JTj), we should note that grain growth is 
assumed to obey a power law (Atkinson 1988, Martin et al. 1997) according to 

R{tr-R'^ = Kt, (11) 

where R is the mean grain radius and K is a constant with Arrhenius temperature 
dependence. In the limit where R{t) ^ Rq, we write R{t) = {Kt)°', with a = 1/m 
being the grain growth exponent. In mean field (non-topological) theories a = 1/2 for 
both d = 2 and d = 3, i.e. a parabolic growth law. Grain growth experiments usually 
exhibit a < 1/2, the difference could be a result of neglecting the topological constraints 
and the detailed microstructure. 



3. General solution 

Following the method outlined in (Forsberg & Massih 2007), we write the general 
solution of equation ([7]) in terms of the basis components 

oo 

C{p,T) = Y,(^i(^^iP,r), (12) 

i=0 
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where are real numbers (constants). Next, we separate the principal time-dependence 
of ej(p, r) by making an ansatz 

ei(p, r) = exp {-Uir) ei{p, r), (13) 

where cUj is a constant and is weakly time-dependent obeying 
Oe ' 

= A'pei + UiCi (14) 



dr 



d' 



2 



d-l Rr 1 d 



(15) 



with A' = — — + \ p ^ . 

P dp^ V p r'^\ dp' 

which is subjected to the boundary conditions: 9ej(0, r)/dp = and ej(l, r) = 0. 
We now express the function Cj in terms of an infinite series 

oo 

ei = ^i\kkei, (16) 

k=0 

where jAfc; i > 1, with jAq = 1, and Ui are selected such that the imposed boundary 
conditions are satisfied for all r and k^i satisfy equation (fT4l) with solutions expressed 
in binomial series in the form: 

.e.(p,r) = ^(^V2r)'^-r/(p), (17) 

where ™/ satisfy the following differential equations 

A;°/ + ^.°/ =0; (18) 
A;r/ + ^.r/ = 2mr'/; form>l (19) 

with ™/(l) = 0, ["/(O) = 5rao and ™/^(0) = 0; Srno the Kronecker delta and fp = df/dp. 
As will be shown below, °/ are expressible in terms of the confluent hypergeometric 
functions, whereas ™/ with m > 1 are related to such functions. 
Writing equation ( fTSl) in a spherical grain, we obtain 

c/^f ,d — 1 , (if ^ „ 

subject to f(l) = 0, f(0) = 1 and fp(0) = 0, where f = ™/, c = dlnTZ/dr = 7^^/7^ 
and w = cUj. Equation fl2Ul) can be transformed to the more familiar confluent 
hypergeometric equation 

d'^u ,d ,du ^ 

subject to u(— c/2) = 0, u(0) = 1 and that Uz{0) exists, where z = — cp^/2, a = w/2c 
and u(z) = f (p). We should note that equations fl20|) and fl2T]) are tacitly time-dependent 
through c. When TZr/TZ =constant, c becomes a time-independent constant and exact 
solutions can be derived. The assumption that c= constant is exact for a parabolic 
growth law, since by recalling equation ([H]), it implies dR'^/dt = 2Dc. 
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A solution of equation (12 ip is the confluent hypergeometric function (Dennery & 
Krzywicki 1995) written in the form 

u(z) = $(a,^;z), (22) 

where ^{a,b]z) is the Kummer function; also with alternative notations M{a,b,z) or 
iFi(a, b, z), see e.g., (Abramowitz & Stegun 1964). This function is regular at z = and 
can be expressed as a power series: 

^ T{a + s)T{b)z' 
^a,b-z) = \ ^ V ^3 

j^o r(a)r(6 + s)r(i + s) 

4. Special solutions 

A point worth noting is that equation 0211) can be expressed as an eigenvalue problem 
of the form Cu{z) = anu{z), where a„ is the energy eigenvalue of the operator 
C = zd'^/dz^ + {d/2 — z)d/dz. We "quantize" the eigenvalues: a„ = + 1/2, where 
n = 0, 1, 2, . . . (n G N) and thus write 

„M(^) = $(n+^,^;^), (24) 
or alternatively 

'j{y) = ^{n + \,^--y% (25) 

Here the choice of a„ = n + 1/2 is merely to obtain a series of convenient expressions 
(harmonics) for the solutions °/(?/) and integrals of lf{y), see below. For the ground 
state, n = and = 3, we have 



'Jiy) = ^ erf(y), (26) 

where y = (c/2)^/^p and erf(?/) is the usual error function. 

Next, returning to equation (1191) and writing it in the spherical system, 

r/pp + + cp) Tf, + ^^Tf = 2m r"V; for m > 1 (27) 

Again, replacing z = — cp^/2, equation (1271) is transformed to 

^X.+ f^ + i-^)X-f^ + i)> = --rV form>l (28) 



n zz ' y ' 2 J z V'^ ' 2/ c 
where the aforementioned quantized energy eigenvalues, ci;/2c = n + l/2, were used, and 
in addition, the spatial dimension was quantized, in the manner, d = 2m + 1. Let us 
write (|28l) for m = 1 and ra G N with the convention „u = ^u, namely 

z nUzz + ~ ^) ~ + ^) " (29) 
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In writing the right-hand side of equation fl29l) . we have utihzed Kummer's relation 

+ i m+^;^) = e^>(^), (30) 

y{z) = <l>{m-n,m + ^;-z), (31) 

where ^{z) is a polynomial when n > m. In particular, with = Ip, have 

ip{z) = l, (32) 
2z 

,p(z) = l + —, (33) 

3K^) = l + y + — , (34) 

4^2 , , 

4p(2;) = l + 22 + — + — , (35) 

Let us now remove from equation (!29l) by making the substitution „m(2) = e^„w(2;), 
we obtain 

Z nVzz + (| + ^) n^'^ - ('^ - 1) n^' = (36) 

This can be differentiated n — 1 times to give 

+ + i + ^-9;*^. (37) 

where the superscript in the parenthesis on nV denotes the number of differentiations 
with respect to the argument. We find „u by first integrating fl37j) once, which yields 



then with further integrations, n times, we obtain 

""l^) = -W^c i - I (7) 

After some rearrangement, we write 

on— 1 rl rz 

"^(^) = -7^^ TTTT / (1 - ^r-'^'e-^' / s"-ie'^^ cisrf«:. (40) 

[Zn — ij!!c Jo Jo 

The integral over the variable z can be expressed in terms of gamma functions, viz. 

s^'e"' ds = (-K)~"-^[r(n + 1, -kz) - nT{n)], (41) 



where r(a, x) is an incomplete gamma function (Abramowitz & Stegun 1964). 
Substituting this result in equation ( HOl) and making use of the properties of gamma 
functions (Abramowitz & Stegun 1964), we finally write 



^viz 



(2n- l)!!c 



Jo ^ ^ ^ ^ ^r(n + j + l) ^ ^ 



i=o 

Hence solution to equation (12^ is found through „m(z) = e^„w(z) 
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Let us evaluate the solution of equation (l^Tj) . relation fl22p . when c is small, i.e., 
TZr/TZ ^ 1 or a ^ 1. We first write fl22|) in a quantized form 

$(a, ^$(a + m, ^ + m;z), (43) 

for m G N. For negative z, we can expand ( |43|) in series (Abramowitz & Stegun 1964) 
according to 

J , oo 

+ ^, I + ^) = r(m + |)e^/2 (I) ' 5^ A„ (i)" J_,(a), (44) 

n=0 

where a = \/ {d — 2m)z — Aaz, u = 1 — m — d/2, Aq = 1, Ai = 0, A2 = d/A + m/2, and 

(n + l)A„+i = {n + d/2 + m- l)A„„i + (2a + m - (i/2)A„_2- (45) 

In equation (jS]), Jp(cr) is the Bessel function of the first kind, related to the spherical 
Bessel function via jp(cr) = a/ vr/2cr 7^+1/2(0"). Specializing to d = 3, equation (jUj) can 
be expressed as 

00 

$(a + m, - + m; ^) = (2m + 1)!! e^/^ ^ A„^V-("^+")j^+„(a), (46) 

n=0 

Since = 0, the first term in series 0461) offers a good approximation for large a. For 
m = 0, we find 

3 , sin a 

2 0" 

with a = ^^(3 — Aa)z. Recalling the substitutions a = uj/2c and z = — cp^/2, we write 

o(^.|-cpV2).;r!:^''''f;f^i. (48) 

/ PV ^ ~ 3c/2 

Hence for c = 0, we have $ = sin (py/cu) /(p^/u)] the result which could be obtained 
directly from equation fl22|) . viz. 

lim$(f ,^;-cpV2) =oFi(;^;-a^pV4). (49) 
c^O /C 2 2 

The function has the series expansion oFi{; b; z) = J2T=o(^i^)/^i^~^^))^'' /^^-^ '^^^ 

be expressed in terms of the Bessel functions of the first kind, e.g., o^i(; 3/2; —ujp^/A) = 

Jo(Pv^) = sm{py/u)/{py/uj). 

Let us now be more explicit and write the series solution from equations f|T3|) and 

( fT6l) as follows 

ei(p,r) = (^oei+ jAi + iX2 2ei H ) exp(-cjir), (50) 

where, according to equation (ITTl) . we write 

oe. = °/(p) 

le. = ,V(p) + 2r°/(p) 

2e. = ^/(p)+4ri/(p)+4r2 0/(p) 

etc. (51) 
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The functions satisfy the differential equations (fTS]) - (IT9]) or (I27|) . For c = 0, i.e. when 

{dR^/dt)/D < 1, we have {d = 3) 

■fip) =Mp^i^) (52) 
Tf{p)=P7'''\pp^rUPP^), (53) 
where /li = ^JuTi. Satisfying the boundary condition °/(l) = 0, we obtain 

/i,=j7r j = 0,l,2,.... (54) 

Next, we calculate the the total number of particles per unit volume of the grain 
at "time" r; i.e. equation (ITUi) . expressed as 

oo 

Q{r) = Y,Q.[T), (55) 

where using equation f[T2]) and in three dimensions, we write 

^i(r) = 3 / x'^ei{x,T)dx. (56) 
Jo 

Substituting for ej(x, r) from equation ( l50l) and employing the relations in ( ISTl) gives 

+ .A2 ( ?/(x) + 4r J/(x) + At' °/(x)) + ■ ■ ■ ] e-'^^rfx. (57) 
The terms in equation (IFTI) can be integrated by using the relation: 

"/(x)rfx = /ir^jn,+i(/^i)- (58) 







In a numerical treatment the parameters iXi,iX2 are determined by satisfying the 
boundary conditions at the beginning and the end of each discretized time interval. 
Moreover, the coefficients in equation ( |T2l) are determined from the the spatial initial 
condition C(p, 0) = Cq. Assuming the orthonormality of the solutions, then in general 

/„^C(p,0) e,(p,0)p-^-Mp 
/oe.(p,0)V-^rfp 



5. Evaluation 

Let us now evaluate the time variation of the number of particles contained in the grain 
(survival rate), viz. 



ds d 



— = Kd-^J r''-'C{r,t)dr, (60) 







Diffusion of particles in an expanding sphere with an absorbing boundary 



10 



where = 27r"'/^/r((i/2). Differentiating under the integral and using equations 
and ([1]), we obtain 



dS 



^ = K.DRitr^ ^-^^ (61) 



dt 



dC{r,t) 



dr 



--m 



We can also calculate the time derivative of G{t), expressed in the form 

where over-dot denotes temporal differentiation. Substituting for S from equation fl6Tl) 
and simplifying, we find 

This equation shows the balance between particle diffusion and grain growth. If we 
designate the characteristic frequency of particle diffusion by z/^ = D / and that 
of grain growth by = R/R, we notice that when -C i'g fhe grain growth is 
dominating and particles cannot reach the boundary of the cell and G = Co{Rq/ RY 
with i?o denoting the initial cell size. Using the relation for the grain growth, R = [Kt)°', 
we can write 

G{t) ~ t-'''^ for ug^ud. (64) 

For example, in 3-dimensional spherical grain with a = 1/2, G{t) ~ t~^'^. Thus the total 
number of particles in the grain remains constant with time, i.e., S = V{t)G{t) =const. 

On the other hand when vd ^ fhe diffusion flux is the controlling parameter 
for escape rate. It may be illustrating to evaluate G{t) in this regime for d = 3 when R 
is stationary, with the well-known exact solution (Carslaw & Jaeger 1959) 

m =3^f^)^ for z.G«l 



oo 



R\ dr J R 

T%2 



G(r)^6C„^ ""— (65) 

k=l 

Here T = Dt/R^ and also equation fl65l) gives G(0) = Cq, as it should. In the long-time 
limit only the slowest decaying eigenmode contributes; hence the asymptotic number of 
particles per unit volume decays according to 

G{T) ~ Coe""'^ for z/^ < 1 A T > 1. (66) 

The corresponding formula for particle survival in the grain is 

S{t) ~ (i^t)3"e-'^'^^""^*'"'" for < 1 A Dt/R^ > 1, (67) 

where again we employed R = (Kt)". Hence for a < 2, S(t) exponentially decays to 
zero at long times. 

Let us finally investigate the case of a slowly growing grain by means of the adiabatic 
approximation (Krapivsky & Redner 1996). In this method, the asymptotic solution 
is constructed to have the same functional form as in the stationary grain, but in 
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addition to satisfy the time-dependent boundary condition, namely, C{R(t),t) = 
and C(r, 0) = Cq. We write {d = 3) 



fit)=Cad{r,t), 



where the function f(t) is determined by substituting equation 
equation ([1]), which yields 



(6^ 



into the diffusion 



fit) 



TT^D + ARR 
4R2 



1 + 



27TrR 



— tan 



/ 7rr \ 
\2r) 



fit). 



(69) 



TT^D + ARR 

Now by placing R = {Kt)°', with a < 1/2, into equation fl69l) the second term in the 
square bracket is negligible and we write 

df ^ ,Tx^D R^ 
Itt 



15)f^ 



4/?2 ' R 

If we further consider that ud ^ uq, then f{t) is given by 

TT 



(70) 



f{t) ^ exp 



-t 



l~2a 



(71) 



4(1 - 2a)K^'' 

which is the same result obtained by Krapivsky & Redner (1996) for the one-dimensional 
case. The corresponding relation for S{t) is 

which is consistent with equation (l67|l . 



-t 



l-2a 



(72) 



6. Summary 

In this note we have treated the problem of particles undergoing Brownian motion 
in an expanding sphere whose surface absorb the particles (Dirichlet condition). The 
time-dependent diffusion equation for particles in dimensional sphere, which accounts 
explicitly for its expansion is formulated. The equation is solved by separating the 
space and time variables and the spatial equation for the particle density is reduced 
to an ordinary confluent hypergeometric differential equation, which implicitly contains 
the grain growth process. For a parabolic grain growth law analytical solutions in 
terms of special functions are derived. When grain growth is dominant over particle 
diffusion, the number of particles per unit volume contained in the grain is calculated 
to evolve with time as G{t) ~ t"""^, where a is the growth exponent and d the spatial 
dimensionality. For the case of slow growth, i.e. when the ratio of grain growth rate 
to diffusivity is small an asymptotic solution is attained. In this regime the long-time 
limit gives G{t) ~ exp[— 7r^(D/A'^")t^~^"], where D is the particle diffusivity and K the 
growth constant. The paper has focused on the mathematical and physical aspects of 
the problem, numerical evaluations for a similar kind of problem are treated elsewhere 
(Forsberg & Massih 2007). 
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